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Abstract 
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an approximate time reversal of the Dyson Brownian motion combined with the convergence 
of the eigenvalue density to the Wigner semicircle law on short scales. 
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1 Introduction 



The fundamental reason why random matrices have been used to model many large systems is 
based on the belief that their local eigenvalue statistics are universal. This is generally referred 
to as the universality of random matrices. It is well-known that the local behavior of eigenvalues 
near the spectral edge and in the bulk are governed by the Tracy- Widom law and by the Dyson 
sine kernel, respectively. Since the seminal work of Dyson [7] for the Gaussian Unitary Ensemble 
(GUE), the universality both for the edge and the bulk were proven for very general classes of 
unitary invariant ensembles in the past two decades (see, e.g. [15l [TBI IH [2l O [6l [M] and references 
therein). For non- unitary ensembles, the most natural examples are the Wigner matrix ensembles 
[19j . i.e., random matrices with independent identically distributed entries. The edge universality 
for these ensembles was proved by Soshnikov [17j using the moment method; the bulk universality 
remained unknown due to a lack of method to analyze local spectral properties of large matrices 
inside the spectrum. For ensembles of the form 

H + aV, (1.1) 

where H is a Wigner matrix, V is an independent standard GUE matrix and a is a positive constant 
of order one (independent of N), the bulk universality was proved by Johansson [13j . (Strictly 
speaking, the range of the parameter a in |13j depends on the energy E. This restriction was later 
removed by Ben Arous and Peche [T], who also extended this approach to Wishart ensembles). 

The approach of [13j is partly based on the asymptotic analysis of an explicit formula by Brezin- 
Hikami for the correlation functions of the eigenvalues of _ff + aV. This matrix can also be 
generated by a stochastic flow ^ 

H + Viy, s > 0, 

and the evolution of the eigenvalues is given by the Dyson Brownian motion |8]. The result 
of [131 H] thus states that the bulk universality holds for times of order one. The eigenvalue 
distribution of GUE is in fact the invariant measure of Dyson Brownian motion. (Rigorously 
speaking, the Brownian motion has to be replaced by an Ornstein-Uhlenbeck process, but we will 
neglect this subtlety.) It is thus tempting to derive the universality of H + ^/sV via the convergence 
to equilibrium. We have recently carried out this approach |12j and the key observation is that the 
sine kernel, as a property of local statistics, depends almost exclusively on the convergence to local 
equilibrium. With this method we have reduced the necessary time to N~^~^^, for any <^ > 1/4 in 
[12] . Note that the relaxation time to local equilibrium is N~^; the additional exponent ^ is due 
to technical reasons. 

From the stochastic calculus, one can see that the typical distance between the corresponding 
eigenvalues oi H + yGV and H is of order {s/Ny^^. Thus the bulk universality of H would hold 
if we could prove the Dyson sine kernel for time s ^ 1/N. On the other hand, for time smaller 
than 1/A^, the eigenvalues do not move in the scale 1/iV and the dynamical consideration seems 
to be pointless. In this gaper, we provide an approach to address the comparison of eigenvalues 
between H + and H. To describe the idea, we now introduce the notations. 

Fix N G N and we consider a Hermitian matrix ensemble oi N x N matrices H = (hek) with 
the normalization 

hek = N^^^^Zik, Zik = xik + iytk, (1-2) 

where x^^, y^/c for £ < k are independent, identically distributed random variables with distribution 
1/ that has zero expectation and variance i. The diagonal elements are real, i.e. ygi = and xu are 
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also i.i.d. with distribution v that has zero expectation and variance one. The diagonal elements 
are independent from the off-diagonal ones. 

Suppose the real and imaginary parts of the offdiagonal matrix elements evolve according to 
the Ornstein-Uhlenbeck (OU) process 

1 X d 

0,ut = Lut, L = -—--- (1.3) 

with the reversible measure ii{dx) ~ dx and initial distribution uq ~ u (strictly speaking, a 
differently normalized OU process is used for the diagonal elements but we omit this detail here). 
Under this process, the matrix evolves as 

t_e-*/2i7+(l-e-*)i/V 

and the expectation and variance of the matrix entries remain constant. Notice for time t small, 
t !v when compared with (|l.ip . after a trivial rcscaling. 

The initial distribution of all the matrix elements is Fd/i*^" — {u d^)*^" with n = N'^ . Let C 
be the generator on the product space and e*'^ := (e*-^)®" be the dynamics of the OU process for 
all the matrix elements. The joint probability distribution of the matrix elements at time t is then 
given by 

i^tdAi®" := e*^u^" d^i®*" = (e*-^u)®" d^i®". 

Suppose that for some t small, say, t = N^^^^ with A > 0, we know the local eigenvalue correlation 
function w.r.t. Ft. Let 



Var{F,Ft) - j \F ~ F^ii®"^ 



be the total variation norm between Ft to F. In order to approximate the correlation functions of F 
by Ft in a weak sense (tested against bounded observables), we need Var{F, Ft) 0. Heuristically, 
Var(F,Ft) ~ tN'^ and this requires that t ^ N^"^ which is far from the time scale t > N^^^^ for 
which the sine kernel has been proven in [T^]. For observables on short scales, an effective speed 
of convergence for the total variation is needed. For example, to test a local observable with two 
variables in scale 1/A^, as in the case of the Dyson sine kernel, one has to prove Var{F, Ft) = 
o{N-^). 

Although the heuristic bound Var{F, Ft) ^ tN"^ can be improved to Var{F,Ft) ~ <iV, fur- 
ther improvement seems to be impossible. Thus we are unable to obtain even the weaker bound 
Var{F^ Ft) = o(l) for t > 1/N. The main observation in the current paper is that, while we cannot 
compare F with Ft, it suffices to prove the existence of some function G for which the correlation 
functions with respect to e*^G can be computed for t > N^^^^ and Var{F,e^''G) = o{N~^). 
Since the necessary input to compute the correlation functions is the validity of the semicircle law 
on short scales, which we have proved for a wide class of distributions v in [SKTUIIII], the choice 
of G is essentially dominated by the condition Var{F, e^'-'G) = o{N~'^). Note that G itself may 
depend on t. Since F = e*^(e~*^F), we could, in principle, choose G = e~*^F = [e~*^]®"F. But 
the diffusive dynamics cannot be reversed besides a very special class of initial data G. However, 
we only have to approximately reverse the dynamics and the choice Gf = [l — tL + ^t^F^] '^"f 
turns out to be sufficient. In this case, e^'-'Gt — F = 0{N^t^) and we will show that 

\Var{e'^Gt, F)\' < [ '"^^^l ~ d^.^" = 0{t'N'). (1.4) 
' ' J e^'-Gt 
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Furthermore, under some mild regularity condition on F, Gt is in the class for which we can 
establish the local semicircle law . We will call this argument the method of time reversal. 

We now summarize the assumptions on the initial distribution. Let the probability measure of 
the real and imaginary parts of the off-diagonal matrix elements be of the form 

i^{dx) = e^'^^^Mx u{x)fi{dx) = e"^(^)e-^'da; 

with the real hmction V{x) ~ U{x) — x^ and similarly for the diagonal elements i7(dx) = e^^'^'da;, 
V{x) = U{x) — ^x^. Suppose that V e C*'(M) and the derivatives satisfy 

6 

"^\V^''>{x)\<C{l + x^)'' (1.5) 

for some fc G N and 

i^{x) < C'e-'^l^l' (1.6) 

with some constants d > 0, C and C". In Section [5] we explain how to relax this latter condition 
to exponential decay, 

iy{x) < C'e-^l"l (1.7) 

with some constants C, C (in fact, some high power law decay is sufficient). We assume that the 
first moment of ly is zero and the variance is ^ 

xdv{x) =0, J x^du{x) = i. (1.8) 

We assume the conditions (jl.Sp . (|1.6p and (|1.8p for V as well with the variance changed to 1. 
Let pn{xi, X2, . • . , xjv) denote the probability density of eigenvalues and for any /c = 1, 2, . . . , TV, 

let 

p^^\xi,X2,.-.Xk):^ j PN{xi,X2,.-.,XN)dXk+l---<^XN (1.9) 

be the fc-point correlation function. With our choice of the variance of the density p^^^ ix) is 
supported in [—2, 2] +0(1) and in the N ^ 00 limit it converges to the Wigner semicircle law given 
by the density 

Qscix) = i-V(4-x2)+ . (1.10) 

Theorem 1.1 Let the probability measure of the matrix elements satisfy conditions (|1.5p . (jl.6p 
and p.Sp . Then for any u with |u| < 2 and for any compactly supported and bounded observable 
O e Lf{M?) we have 



lim / 0(a,/3)7 ]-^P'n{^+ m ^ t ^ + m ^ i N ^dad/g 



0{a,l3) 



1 



sin7r(a-/3)^2n ^^'^^^ 



7r(a - /3) 



dad/3. 
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Remark 1.1 With similar methods we can also prove that the higher order rescaled correlation 
fimctions, 

1 (fe)/ ai 02 Ofc \ 

[Psc{u)]^ r+p,,(u)iV'"+p,,(u)iV--'''+p,,(u)7vJ' 

converge in the weak sense to det(/(ai — where /(r) = ^'"^'^ , however this statement 

requires more regularity conditions on V . The proof of the sine kernel for e^^Gt immediately implies 
the convergence of the higher order correlation functions with respect to the evolved measure. To 
conclude for the higher order correlation functions with respect to F, however, one needs to improve 
the accuracy in ()1.4|1 . This can be achieved by approximating the backward evolution e^*'^ to a 
higher order. For example, using Gt — [l~tL + ^(— tL)^ — . . . (^^l^)! {~tL)"^~^'\ "F, will improve 
the bound (|1.4p to modulo corrections, if V is 2m-times diffcrentiable with bounds 

similar to (|1.5p . 

Remark 1.2 With the same method, the condition that V £ C^(M) in Theorem ll.ll can be relaxed 
toVe C^+^iR) for any e > 0. Hcuristically, this can be seen by observing that, with F ~ w**" 
and Gt = ff " = [l-tL + ^t^L^fF, one can estimate the difference \e*^Vt -v\< 0(1^+" L^+^v), 
which, compared with the estimate 0{t^L^v) used in p.4p gives less decay (still enough to deliver 
the result of Theorem II. but requires less regularity of v (only 4 + 2e derivatives). A rigorous 
proof of this fact can be obtained by using part of the evolution eJ'' to regularize, on the scale t, 
the initial density. 

We now state our result concerning the eigenvalue gap distribution. For any s > and |u| < 2 
we define the density of eigenvalue pairs with distance less than s/N gsdu) in the vicinity of u by 

A(«; x) = ^ #1 1 < J < - 1 : - X, < -u\< tA (1.12) 

2NtNgsc{u) L Nqsc[u) i 

where tpf = N^^^^ for some < (5 < 1. 

Theorem 1.2 Suppose the probability measure of the matrix elements satisfies conditions (|1.5p . 
(|1.6p and (jl.Sp . Let /Cq be the operator acting on L^((0, a)) with kernel ■ Then for any 

u with \u\ < 2 and for any s > we have 

r d2 

lim EA(u;s,x)~ / p(a) da, p( a) = -—r dct(l ~ K-a) , (1-13) 
N^oo Jq da^ 

where det denotes the Fredholm determinant of the compact operator 1 — ICa ■ 

As a corollary of Theorem 1 1.2) one can easily show that the probability to find no eigenvalue in the 
interval [u,u + a/{gsc{uo)N)], after averaging in an interval of size A^^-'^+'' around uq e (—2,2), 
is given by det(l — /Cq,), same as in the case of GUE (see, e.g., [3]). Note that assuming more 
regularity on the exponent of the density u(x) ~ e'~^^^\ we can get a better bound on the 
convergence rate (by approximating the backwards evolution e~*'^ to a higher order) and avoid 
therefore the averaging over u. 

The proof of Theorem 11.11 and 11.21 consists of two main parts. In Section [2] we prove the 
approximation (|1.4p under precise conditions on the initial distribution u = e~^ . In Section [3] we 
prove the sine kernel for the distribution e^^Gt with t = N^"^^^ for any A > 0, which is the optimal 
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time scale lor such a result. Our approach is to recast the formula for the correlation function in 
[T3| . which becomes unstable lor t 1, into a more symmetric form (Proposition 13. 2p so that it 
is stable lor all time up to t = iV^^+'^. The saddle point analysis can then be achieved with the 
local semicircle law from [11]. Finally, we complete the proofs of the main theorems in Section [H 

The method of time reversal described previously is very general and should be applicable 
to a wide range of models. More significantly, it explains the origin of the universality, i.e., 
the universality comes from the "time reversal" . To summarize, the universality consists ol the 
following observations: (1) The local statistics are determined by the local equilibrium measures. 
(2) The relaxation to local equilibria takes place in a short time. (3) The original distribution 
can be well-approximated by the distribution of the Dyson Brownian motion for a short time with 
initial data given by an approximate inverse flow. To implement this scheme, a key input is to 
estimate the fluctuations of the empirical density of eigenvalues in short scales. 

Shortly after this manuscript appeared on the arXiv, we learned that our main result was also 
obtained by Tao and Vu in [18| under essentially no regularity conditions on the initial distribution 
provided the third moment of v vanishes. Some partial results for the Gaussian orthogonal 
ensembles are also obtained and we refer the reader to the preprint for more details. 

Conventions. We will use the letters C and c to denote general constants whose precise values 
are irrelevant and they may change from line to line. These constants may depend on the constants 
in (iLSll-dlll). 

2 Method of Time Reversal 

Recall the Ornstein-Uhlenbeck process from ()1.3p with the reversible measure fi{dx) = ^{x)dx = 
e~^ Ax. Let w be a positive density with respect to /i, i.e. J udfi = 1 and we write u(x) = 
exp{—V{x)). 

Proposition 2.1 Let V satisfy the conditions (|1.5p . (|1.6p with some k and (|1.8p . Let \ > be 

sufficiently small and t = N~^~^'^. Define a cutoff initial density as 



where is a smooth cutoff function satisfying 9{x) ~ 1 for \x\ < 1 and 9{x) = for \x\ > 2 and 
cpf and djv '""e chosen such that v{x)dfj,{x) is a probability density with zero expectation. Denote 
C = i®", F = u'^" and = z;®" with n = . 



v{x) : 




Vc{x) V{x)9{{x - C7v)iV-^/4'^) + d^, 



(i) We have 




(2.1) 



with some c > depending on k and A. 

(ii) gt {1 — tL + ^t^L'^)v is a probability measure with respect to d/i and for Gt ■= [gt 
have 



we 




(2.2) 



where C depends on X and on the constants in (jl.Sp . (jl.6p . 
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In the formulation of this proposition we have not taken into account that in our apphcation the 
diagonal elements of the matrix evolve under a differently normalized OU process with generator 
L = _ with invariant measure /"^dx. This modification is only notational and does 
not affect the validity of the estimates (|2.ip and (|2.2p . 

Proof. From condition (|1.6p the estimate (|2.ip follows directly by noting that the constants cjv 
and (In are subexponentially small in N . For the proof of (|2.2p . we first control the evolution of 
each matrix element under the OU process (|1.3[) . We assume that for the initial density v 



Lv{x) < Aiv{x), L^v{x)>-A2v{x), \L^v{x)\ < Azv{x) 



(2.3) 



hold with some constants positive Ai, A2 and A3. Set gt ~ {1 — tL + ^t^L'^)v for some t > and 
note that gt is a probability density with respect to /i if 



r 

tAi + -A2 < 1. 



(2.4) 



Define 



then 



Vt 



e'^^gt 



l-tL + -t^L^ ]v, 



dtvt 



1 



Note that by the monotonicity preserving property of the Ornstein-Uhlenbeck kernel and by (|2.3I 
we have 

e'^L'^v < Aae'^v < Ase'^^'v, s > 0. 



Here we used the fact that e'^^v < e'^^'-v under the first condition in 
integrating the inequality 



ds 



In particular 



Vf ^ V -\ — 
2 



s'^L^e'^v ds > v( 1 - -t^Age*^' 



assuming p.4p and 
Then 



1 



(2.5) 

which follows from 



(2.6) 
(2.7) 



Vt 



d/x 



10 



- 10 ^ 



(2.8) 



where we used ([2J|) . (|231) and finally ([23 
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Now we consider the evolution of the product density — w®", note that J d/i*^" — 1. 
Applying the same procedure to each variable, we have 

j ^""'^^J'^^"^^ d^i^" < e^^3*«" _ 1 < cAlt^n (2.9) 

as long as A\t^n is bounded. In our application n = 7V^, thus (|2.9p will imply (|2.2p provided that 

A3 < Ct-^/^ (2.10) 



which will also guarantee (|2.7p . It is straightforward to check that the density v{x) satisfies p.3p 
with constants Aj subject to (|2.4p and (|2.10p . This completes the proof. q 

3 Sine kernel for the time evolved measure 

We use the contour integral representation for the correlation functions of the eigenvalues of a 
matrix of the form H = H + aV, where is a GUE matrix [31[T3|. We will apply this result for 
the matrix 

e*^Gt = e-*/2[G, + (e*_i)i/V] (3.1) 

where, apart from a trivial prefactor e~*/^, Gt plays the role of H and a = (e* — 1)^/^ « t^/^. In 
order to be able to use the formula given in Proposition 1.1 of [l3j to analyze H = H + aV, we 
rescale the variance of di' from ^ to ^ + which changes the semicircle law for H = H + aV to 



,iu) -^^^V(l + 4a-.^),. (3.2) 



In particular, the support changes from [—2,2] to [— \/l + Aa^, y/l + 4a^]. Since eventually a will 
go to zero, the condition |u| < 2 in Theorem (jl.ip to be away from the spectral edge changes to 
the condition |u| < 1 which wc assume in the sequel. The semicircle law for H will also change 
from the one given in (|1.10p to 

g^^iv) := -7(13^. (3.3) 
In the rest of this Section we will use (13.311. The main result of this section is 



Proposition 3.1 Let p"^^ he the m-point eigenvalue correlation function for the ensemble H + aV 
defined above and let O : W be a compactly supported bounded observable function. Then for 

any \u\ < 1 and a ]\[^'^/^+^/^ me have 



lim 



f ^, , , /sin7r(ai — Q!,)\ ™ 
/ 0(ai,...,a„Jdet ^] dai . . . da 

Jr^ V Tr[ai-aj) / i,j=i 



{a^-aj) / i,j=i 

Proof. Using Proposition 1.1 of [l3], the (symmetrized) distribution of the eigenvalues x = 
(xi , . . . ,xn) of H = H + aV for any fixed H is given by 
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where y = (yi, . . -yN) is the eigenvalues of the Wigner matrix H with the choice oi S — a^/N. 
Note that 

^ ^ (3.6) 

= E/ ^ 0(iV£i(ii)(a-j^ - m), . . . , A''gi(w)(xi,„ - u)j9s(x,?/)da;i . . .dccAT, 

ii,i2,...,i™ = l 

where E denotes the expectation is w.r.t. the H ensemble. Since O is bounded and the sum 
contains N"^ terms, we thus need to compute the limit of the correlation functions of qs{x,y) in 
the X = (xi, . . . ,xn) variables for a large set 3^jv C of fixed y = {yi, . . . , y^) so that 

P(2/(H)^3;jv)-o(7V-™). 

where y{H) ~ {yi{H), . . . ,?/jv(iJ)) arc the eigenvalues of the Wigner matrix H. We will choose 
yN to be the event that the points y = (yi, . . . , yjv) follow the semicircle law p.3p . The limit of 
the correlation functions of qsix, y) will be computed starting from the next section in Proposition 



More precisely, let 

ri - V? (3.7) 

with some sufficiently small t^o < 1 and we set 



3;.:={,eR-: sup 11^:-^- /^^|<A.-/4 supb|<i^} 

^ Imz>7] ' -^^ J ^ yj J ' J -* 



(3.8) 



for some large constant K. 

By Theorem 3.1 of [11] we then have 

HyiH) ^ Vn) < (3.9) 

(after taking the supremum over all energies, which can be controlled taking energies on a grid of 
spacing r]) . Note that the variance of the matrix elements in [11] was different (see remark at the 
beginning of Section [O]) but this does not change the estimates. The condition CI) of [Tl] on the 
Gaussian decay for the initial density gtfi = {I — tL + ^PL^)vij, is clearly satisfied by (|2.3p and 
(|1.6p . Combining the estimate p.9p with Proposition 13.31 and with the argument after p.6p , we 
have proved Proposition 13. II q 

3.1 Contour integral representation of the correlation function 

We compute the correlation functions of qs{x; y) in a;, for any fixed y G yn'- 

p-j^l c;{xi,...,Xm) = I qs{xi,...,XN]y)dxm+i---dxN- (3.10) 

Note that this definition of the correlation functions differs from the definition of i?^ given in [T3] ; 
the relation being 

r>N I . \ ~(™) / \ 

The following representation is based on the formula in [13j , but it is more stable and suitable for 
analysis for very short time. 
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Proposition 3.2 The correlation functions can be represented as 

RZi^i, ...,Xm;y)^ det{lCff{xi,Xj;y))'!^^^^, (3.11) 

where 



N 

w-Vj 



Z — Vj 

(3.12) 

X -^(w -r + Z-U -SY- \i^^-2u^-.^+2u.)/2S^ 

where r S M is arbitrary and 7 = 7+ U 7_ is the union of two lines 7+ : s — *■ —s + iuj and 
7_ : s ^ s — icj (s G for any fixed w > and F is s — > is, s G M. 

We note that T can be shifted to any vertical line since the integrand is an entire function in 
w and has a Gaussian decay as \Im w\ 00. The constants r e R and uj > Q (appearing in the 
definition of the contour 7 in Kf^) can be arbitrary and will be specified later. 

Proof of Provosition [STB . From Eq. (2.18) in [13], we have 

R^{xi,...,Xm;y) = det{Kff{x^,Xj;y))'^_.^^, (3.13) 

with 

N 



(2TnyS Jr^ w-z ^J-^ z - y. 



where 7 is a contour around all the yj, j = 1, . . . ,N, and F^ is the vertical line M. 3 s ^ L + is, 
for a fixed L so large that 7 and F^ do not intersect. Eq. p.l3p remains invariant if we replace 

Kn by 

K;^(u,«) = e'■(^-")/V"^-^^)/2^i^^(^,^;) = -i— / dz [ d«; ^^^^^ e^^" 
for arbitrary r G K. Here we defined 



,2 



N 



Hv{w) := — ~ vw ^ log(w - yj) . 
The change of variables w = {1 — f3)r + (3w' , z = (1 — /3)r + (3z' leads to 

KS {u,v) ~ I dz' [ dw' ^'''"""'^^ ^[H^[[l-P)r+Pu,')-HM1-P)r+P.'))/S 

{2TTiyS Jr^ w-z 

for every (3. Taking the derivative in /3 at /3 = 1, and removing the primes from the new integration 
variables, we find the identity 



(2myS J;y Jr^ S 



{w - r)H',{w) - {z - r)H'^{z) 



w-z 
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Using that H^Xw) = w — v + •S'X^jLi l/(^ ^ Vj)^ find 
{w — r)H[,{w) — [z — r)H'^{z) [w — r){u — v) 



and thus 



1 



w — z 



{u — v) 



w — z 



w — z 



dz I dz^e-^-")/^ ^"'^ ,(g„W-g„(z))/5 



dz / d,i;e'^(^-")/^e(^-('^)-^"(^»/^ 



1 
5 



{w ~ r) + 



The second term on the r.h.s. is just {v — u)-^JCn{u, v). Therefore 



— [{v-u)IC'^{u,v)] 



1 



1 

(2771)25- 



dz / d«;e'-(''-")/^e(^"('")-^"(^»/^ 



1 



1 



dz / du;e'-(''-")/^e(^"('")-^"(^))/^4 



N 



w — r + z — u - 



^ (z - yj){w - yj) 



Integrating back over u, starting from u, we find that 
1 



{v - u)IC^{u, v) 



N 

dw (^e~i^-^)i-"-'^)/S _ Ij g(t«'-2««,-z"+2«z)/2S -Q 



1 



(w — r) 



AT 



w — r + z — u - 



1 (^-yj)(u'-?/j) 



At this point the contours of integration can be modified; since the singularity l/{w ^ z) has been 
removed, they arc now allowed to cross. This completes the proof of the proposition. q 

Proposition 3.3 Let k > 0. For any sequence y = y^^'' G with the choice S = 1S[~^~^-^ we 
have 



1 



Af-^oo Ng{u) 



a P 



Ng{u)' Nq{uY'') 7r(a-/3) 



sin7r(Q! — 13) 



(3.14) 



uniformly for |u| < 1 — k and for a, (3 in a compact set. Moreover, the correlation functions satisfy 



lim 



1 



_~(m) 

~Pn,v,s 



0L\ 



' + , N = det — 

Nq(u)J V 7r(a 



sin7r(aj — aj) 



uniformly for \u\ < \ ^ k and for ai, . . . , am in a, compact set. 



{ai-aj) J i,j=i 



(3.15) 



Proof. The statement in p.isp follows directly from ()3.14|) and (|3.1ip , so it is sufficient to prove 
()3.14p . We will prove (f3?T4l) in the form 



Nq{u) 



TTT 
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for any sequence m*^^-* with — u*| < C/N and for every fixed with < 1 — k. In order 

to get (|3.14p . we take u''^^ = + a/Ng{u^) with = u. 
Set 

g^g{u^), ^ N-^+^. (3.16) 

From (Hmi), we find 



with 



Mz) = - 2uWz) + 1 ^ log(z - y,) (3.18) 



<7Ar(z,w) = — -fw-r + z-u^^)] - — -Vt ^^-7^^ (3.19) 

^ t(w;-r)^ ^ 7V(«;-r) ^ (u;-?/j)(z-?/j) ^ ' 



hN{w) = i (e-^("'-'')/*^ - l) (3.20) 
with p.l6p . We wiU need the identity 



5Ar(z,u;)== fN[z) + — ■ (3-21) 

w — r z — w 

3.2 Saddle points 

For brevity, we will drop the superscript and denote u^^-* by u in the sequel and we fix |u| < 1. 
We first determine the critical points of /^r, i.e. we solve 

f^{z)^r\z~u) + l-Y,^—^0. (3.22) 
A* ^—^ z — i/j 

This is equivalent to finding the zeros of a polynomial of degree iV + 1. There are — 1 real roots 
and two complex roots, called q^, that are complex conjugates of each other 

fNi<lN) = 0- 

We will work with qj^, the analysis of the other saddle is analogous. Clearly |i?e qN\ < K 

for some large K. 
We can define 

/(2)==;t:(^^-2wz)+ / £>,c(2/)log(z-y)dy 
Js. 

and instead of (|3.22p . we can solve 



f'{z) ^r\z- u) + 2{z - _ 1) = 0. (3.23) 
The solutions of this latter equation (for small t) are given by 



± _ (2t + l)u ± 2tWl + U 



u 



2 



l + 4i 



u 



(1 - 2t) ± 2tiyJ\ - u2 + 0(f), (3.24) 
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and thus in particular 

Im{q^) =±0{t). 



We have 



/"('?) = 7 + 2 



and 

r V 1 — u 

where we also used the equation p.23p for . We set q = q^. 
We need to know that 7^ at the (/at saddle. 



It follows from (|3.8p that for y G 3^ we have 

C 



^supJ/(f)(z)-/W(z)| < je^^, (3.25) 



by contour integration. 

We compare q and ^at. We have fr'om p.22p 



and 



= -Fjv((7jv) := ^4 - ^ , ImqN>0 (3.26) 

N ^ qN ~ Vj 



0sc{y)dy 



q = F{q):=u-t ^^^^^ = u-2t{q-^q^) (3.27) 

First we show that for the only solution to ()3.26p with positive imaginary part we have Im q^ > V- 
This is a fixed point argument. 
Define the compact set 

S := |z : \Re z ~ u\ < Ct, -q < Im z < Ctj 
for some large constant C . Since y e 3^, we know that 

Ct 

sup \Fm{z)~F{z)\ < -3^. 

Im z>rj J V 

For z £ E clearly 

Re F{z) = u + 0{t) 

and 

Im F{z) = 2t\/\ - V? + 0{t^) 

thus 

ReFN{z)=u + 0{t), Im Fn{z) ^ 1t\J\ - v? ^ o(t) 
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so Fn{E) c ^. 

Now we compute, for z S 5, 



N 

(here we used p.25p with i = 2 and observed that = tf'^), and 



F'(z) 



-2i 



1 



with F'(z) = 0{t) if z e S. Thus |i^J^(z)| < 1/2 for z e S, so Fn is a contraction on S and thus 
p.26p has a unique solution, which is qn- 
Comparing the two sohitions, we have 

\qN - g| = \FN{qN) - F{q)\ < Sup |F^(z)||gAr -q\ + |Fjv(g) - Fiq)\. 



Since y G 3^, we get 



\FNiq)-Fiq)\<t 



N ^ qN- Vj 



ct 



thus 



Iqn -q\< 



Ct 



(3.28) 



3.3 Evaluating the integrals 

Using Laplace asymptotics, we compute the integrals in (|3.17[) . We choose the horizontal curves 7± 
to pass through the two saddles q^ = a±ib oi f (see p.24p ). i.e. we set w = 6 (see the definition 
of 7^ after p.l2|) ). The vertical line T is shifted to pass through the saddles, i.e. Re T = a. 
Moreover, if necessary, we deform F in a 0(A^^^)-neighborhood of a so that min^ dist(r, yj) > N^"^ 
and dist(r, a^v) > N~'^] this is always possible. 
Wc split the integrals as follows 

^/Civ(", u+^^-y)=A++ +A+- + A-+ + A" 

according to whether Im z and Ivn w are positive or negative, e.g. 

A^^:^n[ ^/,^(^)5^(,,u,)e^(/«(-)-/«W) (3.29) 

J^± 2m Jp± 2711 

where r_|_ = F n {w : Im w > 0} and F_ = F n {w : Im w < 0}. We will work on A^^, the other 
three integrals are treated similarly. 

The main contribution to the integral A'^^ will come from an e-neighborhood in z and w of 
the saddle point qN ~ Qn- "^^^'^ radius e will be chosen such that after a local change of variable / 
and /tv become quadratic near the saddle. We now explain the local change of variable. 

Since /(z) : C ^ C is an analytic function with f'{q) = and f"{q) 7^ for q = there 
exists an invertible analytic map (f) : z ^ 4'{z) in 

D, :={z : |z-g| <e} 
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with (j}{q) = 0, <j)'{q) = ^tf"{q) such that 



f{z) = f{q) + -[<P{z)Y 



with 

Here e must satisfy 



z) = Vtr{q)i^-<l)i^ + 0{z-q)), Z€D,. 

\r{q)\ 



£ < 



2sup^J/'"(z)| 

we also assume that £ < rj. We wiU choose £ = ct with a smah c, depending on u. We have 

/"(g) = ri + o(i), sup|/"'(z)|<c 



(3.30) 
(3.31) 
(3.32) 



(3.33) 



from the expHcit formula jSSSl), so (|3:32l) is satisfied. Note that (j)'{q) = \/t,f"{q) = 1 + 0{t). 
We have a similar change of variables for /at, i.e. 0Ar with the properties that 



and 



with 



</'jv(gw) = 0, 4''n{in) = ^tfj^) = l + 0{t) 

z G -De, TV ^ {z ■■ \z~ qN\<£} 
(pwiz) = ^Jtf'^{qN){z - (7Ar)(l + 0(z - qn)), 



fN{z) = /w(gAr) + ^[0w(2)] 



2 e L'e.AT. 



(3.34) 

(3.35) 
(3.36) 



This holds if 



£ < 



c\fU<lN)\ 

supz,,.j/;^'(^)r 



For y G 3^, we have f'](,{qN) = [l + 0{N~^/'^)] and |/;^'(z)| < Ct-'^N'^'^ by (IX^ and (0321), 
thus we can choose £ — ct for some small constant c < ^/\ — v?. 

Moreover we have |0jv(z)| < C\z — gjvl for — < ct, so by Cauchy formula \(j)'^{z)\ < C and 
|</'Ar(^)| < Ct~^ for — (j| < ct (maybe after reducing c). The same formulas hold for as well. 
We also have 



Wiq)-q^'A'iN)\ < \Vtriq)- v'^/;^(9)| + ytmq) - ^JtfNiqN)\ < CN-^'\ 

where in the first term we used p.25p and in the second we used \f'iq{z)\ < Ct^^. 
From (|3.31[) and (|3.36p we have 

\<P{z) - cI>n{z)\ < \ct>'{q) - c^'{qN)\\z - q\ + |<^'(gAr)|k " Qn] + C\z - q\' < CtN'^/^ 
and then by contour integration 

\^'{z)~^'j,{z)\<CtN-'/' 



(3.37) 



(3.38) 



for any z with — < ct. Therefore the maps ((> and are C^-close within and both of them 
are C^-close to the shift map z ^ z — q. 
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a i; a+ e/2 a+£ 



a+2e 



Figure 1: Integration contours around the saddle qn = qjf 



We first consider the z integration. Recall that Qn = qti ^ + *^jv from p.24p . We fix a 
small positive constant ci <C 1 and we define the domains 

fl := = X + iy : \x — ajv| > £, |y — ^atI < cie/2| 
rj* := |z = a; + : |x — oatI > e/2, |j/ — 6jv| < ci£/2| 

and 

W:—<z~x + iy : |a; — ajv|<2e, |y — 6jv| < ci|a; — aAr| ^ 

where e = ct. Recall that 7+ was the horizontal line going through q = a + ib, the saddle of /. We 
will deform 7+ to 7^ so that it passes through qjy and it matches with 7+ at the points aAr±2£ + i6. 
Within the regime \Re 2; — ajv| < £, we define 7^^ by the requirement that Im 0jv = along 7^. 
Since (f>N{z) is close to the map z — *■ z — g^r by ()3.36|) . clearly 7^ is almost horizontal curve in small 
neighborhood of q^-, so it remains in W until it reaches the vertical lines \Re z — ajv| = £■ In the 
regime e < \Re z — qn] < 2e, we require that 7^^ matches with 7+ at the points qn i2e + ib and 
it remains in the wedge W. In the outside regime, \Re z — ajv| > 2£ we set 7^ — 7"^, in particular 
j^cWun (see Fig. [J). 

Lemma 3.1 We have 

Re[fN{x + iy)- fN{qN)]> ^^{x-af forx + iyeQ (3.39) 

and 

Re[fN{z) - fNiqN)] >0 for zeWand \Re z - a\ < e (3.40) 

Proof. The second statement p.40[) follows from the normal form p.35p and the fact that for 
z G W we have \Im {z — qN)\ ^ Ci\Re{z — qN)\, i-e. Re{z — qN^ — Oj ^^'^ 'pN is close to the map 
z ^ z — qN i'a W, so i?e[(/)Ar(2;)]^ > for z g VF. 
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For the first statement, we assume x > a, the case a; < a is analogous. We get by explicit 
calculation ^ 

Re f'{x + iy) > —{x — a), for x + iy e fi*, x > a. 



2t 



Using (|3.25p for £ = 1, we have 



Re d^fN{x + iy) > Re f'{x + iy) - CN~^/* > ^{x - a) for x + iy e il* x > a (3.41) 

(the error is absorbed since |a; — a| > ct/2 for x + iy £ ft*). Since i?e[/jv(z) — fNilN)] > on the 
vertical lines |a; — a| = e/2, \y — b\ < cie/2, we can integrate the inequahty (|3.4ip to obtain p.39p . 

□ 

In order to estimate the w integration along r+ parametrized as a + is, s > 0, we analyze the 
behaviour of Re f along r+. For |a; ~ a\ < Ct and y G E we first compute 

Re dyfNix + ly) = -Im f'j^{x + ly) = -Im f{x + iy) + 0{N~^^^) 

which holds for \y\ > rj. By explicit computation, and using /'(a + ib) = 0, 

-Im f{x + iy) = -{y-b)[^ + 2^+0{t)+Oiy^) 



if \y\ < ^Vl — u^, \x — a\ < Ct for some large C. Thus we have 



Re dyfNix + iy) < - 



V<\y\<^\/l-u^ and y~b> e/2, 



where e = ct with a small c as before and a similar lower bound holds for y — b < —e/2. Defining 
?1 ■■= = x + iy : e < \y - 6iv|, V < U < ^ Vl - u^, \x - a^l < cie/2^ 

W :^ ^ x + iy : \y - BnI < 2e, \x - unI < ci\y - 5Ar|| 
analogously to W before, we easily obtain 



and 



Re [fN{x + iy) - fN{qN)] < -TZjiv - b)"^ for x + iy e fl 

lot 



Re [fN{w) - fN{qN)] < w eW, and \Imw - b\ < e, 



(3.42) 
(3.43) 



similarly to the proof of Lemma 13.11 

The regimes < y < rj and y > ^\/\ — v? are treated directly. We use 



RedyfN{x + iy) 



Im t \z-u) + ^Y.— 
'J 



(3.44) 



> -v/t. 
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Hence for < y < rj we have 

Re [fNix + ly) - /Ar(gjv)] < Re [/^(.t + z??) - fNiqN)] + |^ ^ "^(^ " 

from p.42p . if 770 is sufficiently small, see (|3.7p . 
liy> i^l -u2, then 



ly L 



< 



N^{x-y,Y + y^ - 713^' 
hence 

Re dyfNix + iy) < -y/2t (3.45) 
and thus Re fN{x + iy) < —y"^ /At in this regime. Summarizing these results, we have 

Re [fN{x + vy) - /jv(gjv)] < -^^(y - hf (3.46) 
holds for any G K and |a; — a| < cie/2. 

We can define a new contour similar to the 7^. It follows the path where 4)n has zero 
imaginary part when \Im w ^ h\ < e/2 and then it returns to F"'" when \Im w — h\> e. We recall 
that miuj dist(F^, > N^"^ and dist(F, ajv) > N^"^ by the choice of F. 

With the paths 7^ and F^ defined, we can now do the integration 

A++:=n[ p-f pLh^{^,)g^^z,n^)e^(M^)-M^)). (3.47) 
•^7+ 7r+ 

Near the saddle we need the bounds 

\gN{z,w)\ < C/t, \d,gN{z,w)\ < C/e, \hN{w)\ < C (3.48) 

if |z — (a + ih)\ < e, \w — {a + ib)\ < e. In order to make sure that these bounds are satisfied, we fix 
the constant r = Re gAr(w*) in (|3.17p . Here qN{u*) is the unique solution with positive imaginary 
part of the saddle point equation (|3.22p . with u (which is actually a short hand notation for u'^-*) 
replaced by the fixed u*. Note that, since \u^^^ — u*| < C/N, we find that the real part of the 
exponent of hniw) (see p.20p ) is bounded, \r — Rew|/ip < C, as w runs through F. 
This choice also guarantees that, away from the saddle, 

|/iw(w)| < C7e^*"'l^""'-"l, \gN{z,w)\ <CN'^ (3.49) 

that hold for |/to z| > 77, /to w > 0. These bound follow from (PTTg)) . ((5?^ and ([X^ and when 
w is near the real axis, we also used that Fjv is away from the j/j 's. 

The integration in (see p.47p ') will be divided into regimes near the saddle gjv ("inside") 
or away from the saddle ( "outside" ) : 

= Aii + AiQ + Aqi + Aqq. (3.50) 
Recall that \qN — q\ = o[t) and q ^ q^ = a + ib (see p.24p ). For example 

A,o-N [ x{Rez-a)^[ (1 - x(/m - 6))l^/i^(H5^(z, «;)e^(^"('")-/"(^». 
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where x is the characteristic function of the interval [— £, e\. The other A's are defined analogously. 
Using ((339)) lOO)) and ([Qe]) . we have 

\A,o\<N I dzx{Rez~a) ( (1 - x(/m w - 6))|5Ar(z, e^-"" [^"('"^^^"('^"^Idw. 

The integral of the exponential term is bounded by 



■Aft 



Taking into account (|3.48p and (|3.49p . we see that \Ai, 



< e~=^* since t = N-^+^. Similarly 



we 

can bound all other terms with an outside part. When \Re z — a\ > ct ^ N^'^, then the exponential 
growth oi hpf in (j3.49p will be controlled by the Gaussian decay of 

g-JVKe[/«(2)-/jv(<2«)] < ^-cNt-^\Rez-a\^ 

from (12311). 

Finally, we have to compute the contribution of the saddle, i.e. the term An. We let 7 be the 
part of 7^ with \Re Jn — o,\ < e and similarly defined F. Recall that Im (j)^ = on 7. From 
standard Laplace asymptotics calculation, we have 



J 'y 

hN{w)gN{qN,w) + n{w) 



2n 



(3.51) 



using (|3.34p with 
Using (|3.48p . we have 



\Q{w)\ < C\l — max \dzgN{z,w)\\hN{w)\ 

1\ zEDc 



C 1 



1171 < ct-\ — 



N t ym 

while the main term in the bracket on the r.h.s. of p.5ip is of order t~^. Analogously performing 
the dw integration, we obtain that 



Ai. 



-1 



'27Tfj(jiqN) 



9N{qN,qN)hNiqN) 



1 + 



-hNiqn) 



2-K 



1 + 



NtJ 



where we also used gN{qN,qN) = f'MiqN) following from (|3.2ip . So far we considered the saddle 
qN ~ qtf with positive imaginary part for both the z and w integrals. The same calculation can 
be performed at the saddle z = w = q^. The mixed case, when z is integrated near one of the 
saddles and w is near the other one, gives zero contribution, since gN{q]^,qjf) = 5iv(9jv, 'Zjv) ^ ^ 
by (|3.2ip . Adding up the contributions of the two relevant saddles, z ~ w — qj^ and z = w ^ q]^, 
taking into account the opposite orientations of the two pieces of 7Ar, one obtains 



1 r 

2^ 



hNiqti) + f^Niq^) 



1 

27rT 



r)/tg , 



r(9iv"'")/*e 



where we used the choice r = Re q%{u*) (see after p.48p '). which guarantees that |r — Re(7^| 
as iV ^ 00, and the equations p.l6p . p.24p . and p.28p . This completes the proof of Proposition 

□ 
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4 Proof of the main theorems 



Proof of Theorem Wc follow the notations of Proposition 12.11 In Proposition 13.11 we have 
shown that the sine kernel holds for the measure e*^Gt ii t = N^^^''^. More precisely, let PN,t{x), 
denote the density function of the eigenvalues x = {xi, . . . ,xn) w.r.t. e^^Gt and let p^^\ be the 
two point correlation function, defined analogously to (|1.9p . Similarly, we define Pn,c{x) and p^v c 
for the eigenvalue density and two point correlation fimction w.r.t. truncated measure F^. = v®" . 
In Proposition 13 . 1 1 we showed that 

lim / \pf(u+^,u+^)0{a,P)daAp^[ 0(a, /3) [l - ( ^^"^^" ~fY j dad/3. (4.1) 

for any \u\ < 2 and with the notation q = Qsc{u). (Wc remark that p\q ^ was denoted by p]y in 
Proposition 13.11 and the condition < 2 is translated into juj < 1 after rescaling.) 
To prove (|l.lip . wc thus only need to control the difference as follows 



(2) 

Pn 



P \ (2) 



where 



{I) 



(2) 

Pn 



NqJ 



Nq'""^ Ng) 



0{a, f3)dad(3 



< (I) + (11), 



a f3 \ 



ill) : 

Using (|2.ip . we have 



(2) / , " I ^ ^ 

Pn'c['^+ — '""I ' 



(2) 

-Pn.c 



(2) 
-PNt 



a (3 \ 

a P \ 

u H , u H 

Nq' NgJ 



Ng' NgJ 

(I) < N^O\\^ J\F- F,|dM^" < Ce-^^' 
with some c > as iV ^ oo. To estimate (//), we have 
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0{a,f3)dadP 
\0{a, f3)\dadp 







(//)</ 



(2) 

Pn., 



(2) 
PN.t 



a 



NgJ 



(2) 

Pn, 



t("+^,«+^)|0(a,/3)|dad/3 



< 



(2) 

Pn.c 

(2) 
PNt 



a P \ , 

u H , u H — 1 

Ng NgJ 



pi^^i(^^+^,"+^)|0(a,/3)|d«d/3" . 
Using (|4.ip for the observable \0\ instead of O, the second factor on the r.h.s. of (|4.2p is bounded. 



1/2 



1/2 



(4.2) 
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Since O is bounded, the first factor is smaller than 



C 



1/2 



< c 



< c 



N 2 

11 pAr,f(a;)da; 



1/2 



Gt - F, 



-I 1/2 



-d/i 



(4.3) 



Here in the first step we used that the quantity D{f,g) ~ J\f/9 ~ Ipff for two probability 
measures / and g decreases when taking marginals. In the second step, we used that D[f,g) 
decreases when passing the probability laws from matrix elements to the induced probability laws 
for the eigenvalues. Finally, we used the estimate (|2.2p . This completes the proof of Theorem ll.il 

□ 

Proof of Theorem \1.2i We first prove Theorem ll.2l for the ensemble H+aV with a = 
(see the beginning of Section [5] for the necessary rescaling). Let E denote the expectation with 
respect to this ensemble and let Ey denote the expectation with respect to the density x — > qs{x,y) 
for any fixed y and S = a'^ /N ~ iV^^+'^. Then we have 

(4.4) 



(4.5) 



E A(u; s, •) = J Ey A{u; s, ■)l{y G y)dP{y) + J Ey A(u; s, ■)l{y e :y")dP(y) 
by recalling (|3.5p . The second term can be estimated by using |A| < N and 



E,y A{u; s, ■)l{y e y^)dPiy) < CNe 



For the first term in (14. 4p . we use the exclusion-inclusion principle to compute 



Ey A(u; s, ■ 



1 X— ^ f f r s 1 

T.yr d„. . . d.„,i{ ».» I.. _ .,1 < -} 

X [jp]y^y^si^ + Vi,U + V2,...,U + Vjn) 



(4.6) 



with g = g(u) (see p.2p ) and recall that p^^y g denote the correlation functions of qs{x,y) (see 
p.lOP ). After a change of variables. 



Ey A(w; s, •) 



1 



2NtNg 



m=2 



Ngt„ 



Ngti, 



Zl , Z„i\. 
11. -i _ 



Np'"' ■ NpJ' 



2Nt 



Ngt 
Ngt 



1 1 max I — I < s I 

-t r^^ g^N 

— V("l)™m / dzi / da2... / da™ 

m=2 J-NgtN Jo Jo 



(4.7) 



1 



mj (NgY 



M - 



Zl + 02 



Zl + a„ 
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where the factor m comes from considering the integration sector zi < Zj, j > 2. Taking N oo 
and using Proposition 13. 3i we get 

^ w ^ (-1)" r, r, , /sin7r(a, -a,)\™ , , 

hm Ey A{u; s, • = 7^ '-ttt / da2 . . . / da„ det , 4.8 

where in the last determinant term we set ai — 0. The interchange of the hmit and the summation 
can be justified by noting that the exclusion-inclusion principle guarantees that (|4.6p is an alter- 
nating series where the difference between the sum and its M-term truncation can be controlled 
by the (M -I- l)-th term for any M. We note that the left hand side of (|4.8p is p{a)da, where 
p{a) is the second derivative of the Fredholm determinant det(l — /Cq.) (see (|1.13p ). Combining 
with the estimate (14.51). we have 



lim EA(u;s,-)= / p{a)da. (4.9) 

After rescaling (j3.ip . we also conclude that the limit of the expectation of A with respect to the 
time evolved ensemble e^'-'Gt (see Proposition [23]) is given by right hand side of (|4.9p . 

Finally, the difference of the expectation of A with respect to the measure e^^Gt and w.r.t. the 
initial ensemble F vanishes since |A| < and Var(e*^G't, F) < CN-'^+^^ (see and (j^ ). 

This completes the proof of Theorem 11.21 q 



5 Some extensions and comments 

In this section we explain how to relax some of the conditions on the initial distribution v. 

We first explain how to extend our proof to include distributions v with compact support. Take 
for example a density w.r.t. the Gaussian measure d^(a;) = that is given by a nice bump 
function u{x) supported in [—1, 1] decaying like (1 ± x)™ near the boundary x = ±1. Clearly, for 
any m fixed this distribution violates the assumptions of Theorem 11.11 We now show that for m 
large enough, it is still possible to prove the universality. Define a new distribution with density 

t"' + u{x) 

li^) = ^^^m (5-1) 

with a small parameter t > to be determined later. Near the edge 1 we have Lq{l — y) < Cy™^^ 
for < y ^ 1 with some m-dependent constant C. We thus need the condition 

C?/"-2 < ^-lyn _^ ym]^ < ?/ < 1, 

to guarantee that (1 — tL)q is a probability density. This inequality holds if 

> Ct. 

The other conditions concerning and (see (|2.3p ) can be handled similarly. Choosing r — 
Ct^/"^ , the total variation norm is bounded by 



22 



Since n — N"^ and t — N ^+^, we have 

Let, say, m > 9, then the error term will be smaller than N^^^^ with some (5 > and this will 
imply Theorem 11.11 for the initial distribution u. The modification of u in (|5.1[) can certainly be 
more sophisticated to reduce the exponent to. 

Second, we show that the Gaussian decay condition (|1.6p can be replaced by the exponential 
decay (|1.7[) . For any ^ > define 



where ai and Z( are chosen so that 



X dvi = 0, / Ave.{x) = 1. 
Due to the assumption (|1.7[) . we have 



\at\ + \Zi -l\ <e 



Let vi{x) ~ £i'e{x£). Clearly, the random variable x distributed according to is bounded by 
1. in particular it has a finite Gaussian moment. Denote the variance of by (t| and we have 
ai = l/£+ 0{e^'^^). We will neglect all the exponential small terms 0{e~'^^) and assume ai = l/£. 
Similar cutoff and rescaling applies to the distribution of the diagonal elements. 

Consider the random matrix generated by the measure vg and i^i and denote the probability 
law of the eigenvalues by fi and fi. Since all quantities introduced below can be defined w.r.t. to 
both and we will only give explicit definitions for !/£. Recall the Stieltjes transform of the 
eigenvalue distribution w.r.t. i^g is defined as 

TO,=TO,(Z)= / (5.2) 

Jr E- z 

where Fi is the empirical distribution function of the eigenvalues. Then the empirical density of 
eigenvalues and rhi converges to the rescaled semicircle law 

pHx) = ipscixl) = ^V4.-xH^, toL(z) = iruscizi). 

We now follow the proof given in [Tl] to prove the local semicircle law Theorem 4.1 [11] for z/^. 
The key estimate is contained in Proposition 4.3 which depends on Proposition 4.5. The random 
variables bj in Proposition 4.5 are now distributed according to and the only assumption of this 
proposition, the Gaussian bound (1.3) (i.e., the condition Cl) is now trivially satisfied since l/i 
has compact support. Hence we can now prove Proposition 4.3 using the same strategy. Thus the 
equation for the probability estimate appearing after (4.6) in the paper still holds but the upper 
bound on the constant defined in (4.6) now becomes 2M£'^ /{Nt]) due to the scaling. Thus the 
key estimate at the end of the proof of Proposition 4.3 of [11] is now changed to 



b[|^| > ^] < 4exp ( - cmm{5^/N^/£, S^Nt]/£^}) . (5.3) 
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Therefore, Theorem 4.1 of [TT] holds with the estimates taking the form 

Pz;i sup \m'^{E + ii^) - mi^{E + tT])\ > s] < Ce-^^^/' (5.4) 

for any 6 < cik/£. Passing from Vi to i^e via scaling, we have 

sup \m^r^{E + ir])-mi^{E + ir])\>5\<Ce~^^V^^ (5.5) 

Ee[-2+K,2-K] 

for 5 < ciK. Comparing this estimate with the original bound (4.1) in (llj . note that the only 
change is that the 77 in the exponent has deterioriated to t]/£. This is due to fact that we applied 
the Proposition 4.5 of [TT] without taking the advantage that the variance is now reduced to l/i'^, 
which should enhance the large deviation estimate (|5.3[) . For our case, however, the estimate (|5.5p 
is already sufficient since we are interested in the case Nrj = iV^ and £ = (log A^)^. 
Finally we need to pass estimates to the original measure vf". We can check that 

Var(t'f",^®") < Cne-"'^. 

Since in our application n = iV^ and £ = (log A^)^, the right hand side is smaller than any negative 
power of N, all necessary expectation values of observables w.r.t. i^f" can thus be passed to ly^". 
This shows that the local semicircle law holds on scales 77 > N~^^^ for any e > assuming only 
exponential bound (|1.7p instead of the Gaussian bound (|1.6p required in CI) of [TT]. This input 
is sufficient to conclude the proof in Section [3] if t = is changed to N"-^'^^ in (I3.16p . 
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